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We have studied the three-band Peierls-Hubbard model describing the Cu-0 lay- 
ers in high-Tc superconductors by using Lanczos diagonahzation and assuming in- 
finite mass for the ions. When the system is doped with one hole, and for A (the 
electron-lattice coupling strength ) greater than a critical value, we found that the 
oxygens around one Cu contract and the hole self-traps forming a lattice and elec- 
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tronic small polaron. The self-trapped hole forms a local singlet analogous to the 
Zhang-Rice singlet in the undeformed lattice. We also studied the single-particle 
spectral function and the optical conductivity. We have found that the spectral 
weight, in general, is similar to that found in previous studies where the coupling 
with the lattice was absent. There is an anomalous transfer of spectral weight but, 
contrary to those studies, the spectral weight goes to these localized polaronic states. 
However this polaronic shift does not seem enough by itself to explain pinning of the 
chemical potential observed in real materials. The peaks in the optical conductiv- 
ity are also shifted, according to the polaronic shift, in the single-particle spectral 
functions. We compare our results to those obtained in inhomogeneous Hartree-Fock 
calculations and we discuss their relation with experiments. 
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I. INTRODUCTION 



Since the discovery of high-T^ superconductors , a considerable effort has been made 
to determine the nature of hole-doped states in the copper-oxide planes which are present 
in these materials. Many theoretical discussions have emphasized the importance of strong 
electronic correlations within these planes, leading to the study of purely electronic two- 
dimensional models like the one-band Hubbard model , the t — J model and three-band 
Hubbard models In spite of the fact that the studies on these purely electronic models 
are still inconclusive, there are clear experimental indications that, at least in the dilute 
regime, the coupling to the lattice should be taken into account p. 

The nature of doping states with and without coupling to the lattice has been studied 
in inhomogeneous Hartree-Fock (HF) calculations for one-band and multi-band Hubbard 
models 0-0. In the three-band Hubbard model, with no coupling to the lattice, a doped 
hole localizes mainly in the four surrounding O sites of a Cu, nearly canceling the spin 
density at the central Cu p|-|r3|. This charge and magnetic polaron is a manifestation, at 
the HF level, of the singlet bound state that Zhang and Rice found in their derivation of 
the t — J model (as an effective Hamiltonian of the strongly interacting three-band Hubbard 
model). By including the coupling with the lattice P,P^ it was found that the O's around 
the central Cu forming the electronic polaron relax with a symmetric pattern [Fig. 1(b)]. 



Similar ideas were proposed earlier in Ref. [jT4|. The infrared (IR) spectra of this state 
PJl3[| , for moderate electron-lattice coupling, calculated in the random phase approximation 
(RPA), shows a doped- induced bleaching of the in-plane Cu-0 stretching mode and a side 
phonon band at lower frequencies. The nature of the charge carriers and their interactions 
with the lattice has been probed by infrared optical absorption in photodoped [p!5|-[T7|| and 
chemically doped samples. Both kind of experiments show in La2Cu04, among other 
materials, a doping-induced bleaching of phonon modes and, in general, an intensity shift 
to lower frequencies. Of particular interest in La2Cu04 is the bleaching of the 708 cm~^ 
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Cu-0 stretching mode and the appearance of a side band at 640 cm~^ |T^JT^. Both are 
clearly resolved in the photodoped spectra and consistent with the data of chemically doped 
samples [|18|. This behavior agrees wells with the HF plus RPA calculations and provides a 
strong indication of the existence of the electronic and lattice polarons. 

In spite of this success of HF plus RPA calculations, the above studies present various 
shortcomings. Since the electronic polaron is already self-trapped without coupling to the 
lattice, a negligible small electron lattice coupling strength A (defined below) suffices to 
produce lattice relaxation. In other words, there is no competition between the delocalization 
energy and the self-trapping effects in this type of calculations. Another serious shortcoming 
comes from the fact that the Zhang-Rice wave function can not be represented by a single 
Slater determinant, and ( at least at the HF level) a relevant part of the physics is missing. 
In this sense, it is important to study the interplay between electronic correlations and 
lattice defects or deformations without the above limitations. The electronic correlations 
can stabilize the lattice deformations as a consequence of the pre-existing enhancement of 
the effective mass. This effect has been observed in an exact diagonalization study of the 



Holstein-Hubbard model [113. On the other hand, the lattice deformations can introduce 
anharmonic effects, further enhancing the effective mass and changing the properties of 
carriers with respect to those found in the framework of purely electronic models. 

In this paper we study, by means of Lanczos diagonalization, the formation and properties 
of polarons in the two dimensional (2D) three-band Peierls-Hubbard model, which we believe 
is a realistic microscopic model of the Cu-0 planes. The model is presented in Section II. 
We take the mass of the O ions equal to infinity (adiabatic approximation) for simplicity. 
This method is free from the above mentioned limitations of HF approximations and takes 
into account all the advantages of exact diagonalization studies of the three-band Hubbard 



models pO|pl| which have reproduced fairly well the main expected spectral signatures. We 
discuss finite size effects in the conclusions (Section VII). By minimizing the total energy, 
we study the range of electron-lattice coupling at which the various lattice deformations 
become stable. 



For any finite electron lattice coupling strength A below a critical value (Ac), we find that 
the system forms a charge density wave (CDW) with a breathing extended lattice distortion 
(BED) [Fig. 1(a)]. A similar distortion was found above a critical value in Monte Carlo 
simulations |22[ of the Hubbard model. In Section III we discuss the phase diagram in the 
parameters space and give general arguments of why a BED should be found for any finite 
A. 

Above Ac a polaron is form. In what follows, we use the term polaron to designate the 
composite object made by the lattice deformation and the self-trapped carrier. The polaron 
results from static oxygen displacements as shown in Fig. 1 (b). This is in accordance with 
the HF picture, but in exact diagonalization calculations, the polaron appears only above 
Ac whereas, as explained before, in HF it appears for any finite A. 

The polaron and the CDW are characterized by studying static magnetic and density- 
density correlation functions (Section IV). The presence of the self-trapped carrier manifests 
also in the single-particle spectral function. The evolution of the states formed by these self- 
trapped carriers as a function of A is study in Section V. Special attention is given to the 
issue of anomalous transfer of spectral weight and evolution of the chemical potential with 
doping pB|J2^ , in view of the controversial experimental situation Finally, the effect of 



the self-trapped states on the optical conductivity is analyzed in Section VI. Consistently 
with the features of the single-particle spectra discussed above, we find a shift of the peaks 
to higher energies. 



II. THE MODEL, PARAMETERS AND CLUSTER STUDIED 



We consider the following Hamiltonian for the Cu-0 planes: 



+i:\ku] (1) 

i 

Here, cl^. creates a hole with spin a at the site i in the Cu dx2_y2 or the O Px,y orbitals, 
^j,(T = c|,o-Cj,(T and rij = (rij^i + nj j). We have included on-site repulsions Ui = Ud,Up for 
holes on Cu and O atoms and a nearest-neighbor Cu-0 repulsion Uij = Upd- Cj = e^, Cp are 
the site energies of the Cu and O atoms. 

The static lattice deformation and the electron-lattice coupling are considered in the 
following way. Each of the patterns of deformation considered here (see next section) are 
defined by a set of displacements Uj of the oxygen ions along the Cu-0 bonds. We assume 
that the nearest-neighbor Cu-0 hoppings are modified by the O- ion displacements as tij = 
tpd ± auj, where the -|- (-) applies if the bond shrinks (stretches) with positive Uj. In order 
to study the stability of the system with respect to a given deformation of the lattice, we 
included the elastic energy in the Hamiltonian (last term of Eq. |I|). For a given value 
of A, and for a set of displacements Uj, we diagonalize exactly the electronic part of the 
Hamiltonian using the Lanczos algorithm. Then the total energy is minimized by varying 
I Uj I constrained to the specific pattern. 

As a reference parameter set we take A = Cp — = 3, t/^ = 8, t/p = 3, Upd = 1, tpd = 1. 
For simplicity the 0-0 hoping tpp is taken equal to zero. We take infinite masses for the O 
ions, so that their kinetic energy drops out. Following the estimates made in Ref. P,|13 



we 



take K = 32tpdA ^. In terms of the above parameters, we have defined the electron-lattice 
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coupling as A = -j^j—- 

The exact diagonalization study was made on the CU4O8 cluster (12-sites) with periodic 
boundary conditions. We consider two cases, the half- filled system (one hole per CU-O2 
cell), and the system doped with one extra hole. 

III. STUDY OF THE STABILITY OF VARIOUS LATTICE DEFORMATIONS 
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As a first step in the analysis of tlie electron-lattice coupling, in this section we study 
the stability of an homogeneous distortion of the lattice (BED) [Fig. |l](a)] and the stability 
of a local contraction of four oxygens toward one Cu (breathing "localized" deformation or 
BLD) [Fig. |l|(b)]. We shall also refer to this localized deformation as lattice polaron. 

In the half-filling condition the lattice remains stable upon these breathing deformations, 
for all physically reasonable values of A (A<7). A similar result was found in a Monte- Carlo 
study of the analogous one-band model [^ and in the HF calculations . 

In general, the system can accept a lattice deformation if it gains enough local covalency 
to compensate the cost of elastic energy, delocalization energy, magnetic energy or Coulomb 
repulsion. At half-filling, the delocalization energy is obviously not relevant. We should 
analyze only the magnetic energy. Coulomb repulsion and covalency. Without deformation, 
the spin density at each Cu site is uniform and equal to a fixed value, say m. When the 
oxygens are displaced, the hopping integrals are modified and it is likely for charges to spend 
more time around the Cu sites of one sublattice and less time in the Cu sites of the other 
sublattice. Thus, for both Cu sites, the spin density m is reduced because, in one case, there 
is an increased tendency to double occupation, and in the other, the sites are becoming 
increasingly empty. The main increase in the energy (of order Ud) comes from the double 
occupancy and a large value of A is necessary to compensate it. More complicated patterns 
of deformation, not studied here, can gain some covalency without increasing the double 
occupancy too much |T2|,|9[|. and, presumably, can destabilize the lattice for smaller A. 

The situation changes upon doping. With one doped hole, the BED is stable for every 
nonzero A smaller than a certain critical value Ac ( Ac ~ 1 for our set of parameters). The fact 
that the BED is stable for arbitrary small A can be understood as follows : the ground state 
of the undeformed system is four-fold degenerated with momentum k = (vr, 0) or A; = (0, it) 
and Sz = ±1/2. Consider now the non-interacting case {Up = Ud = Upd = 0) with the 
same number of particles. The ground state is formed by five-particle Slater determinants 
which, for fixed Sz, can be labeled by the total momentum as, |(7r, 0) > or |(0, tt) >. We can 
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construct two states |+ >= ■^{\{7i,0) > +|(0,7r) >) and |- >= -^^(Kvr, 0) > -|(0,7r) >). 
Taking the x coordinate in the direction of the Cui-Cu2 bond in Fig. |1| it is easy to see that 
1+ > will have more charge in sites 1 and 4 and |— > in sites 2 and 3. Clearly the BED 
distortion favors one of these configurations at the expense of the other. In other words, the 
BED lowers the symmetry of the lattice lifting the degeneracy. This conclusion is valid for 
the interacting case because the ground state in this system has the same symmetry and 
quantum numbers as in the noninteracting one. 

For one doped hole, we show in Fig. |^ the ground state energy, measured with respect 
to the energy of the undistorted lattice, as a function of A. Both localized and extended 
deformations are shown. 

For A > Ac, the BLD becomes stable. The loss in elastic. Coulomb and delocalization 
energies is compensated by the gain in local covalency. In order to understand this crossover, 
we show in Fig. ^, the equilibrium lattice displacements of the oxygens along the Cu-Cu axis 
as a function of A for both types of deformations. In the region A < Ac the displacements of 
the ions are very small and, consequently, there are little effects of the deformation on the 
ground state properties. When A becomes larger than Ac the four displaced oxygens of the 
BLD undergo a sudden contraction toward their central copper. At this point this localized 
deformation becomes stable and the extra hole is trapped in this CU-O4 cluster, forming a 
small polaron state. 

We also studied the variation of Ac when Ud is moved away from the reference parameter 
set (see Fig. p. It can be seen that Ac increases as a function of Ud roughly linearly in the 
studied range {6 < Ud < 20). This is because BLD becomes stable by localizing the added 
hole (see next section) which, in turn, leads to an increasing double occupancy of the central 
Cu site costing an energy of order Ud- Thus, it is necessary a larger electron- lattice coupling 
to overcome this Coulomb repulsion. 

It is worth noting that, in the strong coupling limit, this modification of the hopping 
integrals due to deformation leads also to a local reduction in the effective antiferromagnetic 
coupling J. In fact, in a perturbative calculation, J results proportional to the fourth power 
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of the hopping integral t and then, upon deformation J {t^ — [auYY < f^. The loss in 
magnetic energy, due to the lattice deformation discussed above, could also be related to 
this reduction in J. 

IV. STATIC CORRELATION FUNCTIONS 



The formation of the polaron when one hole is added to the system can be determined 
by analyzing, the spin and charge densities. 

In Fig. 1^ we show the sum of the occupancies < ra, > in the central Cu (Cui) and its 
four surrounding O involved in the BLD, as a function of A. The total occupancy of this 
CU-O4 cluster is ~ 2 for A > 1 indicating the trapping of the extra hole. We also show the 
occupancy of the CU4 site, the most distant Cu from Cui, and its four surrounding O. 

In Fig. ^ we display < Si > m. Cui, CU4, and in the two remaining Cu sites (Cu2 and 
Cus), as a function of A. The reduction of the magnetic moment at the central Cu results 
from the fact that the two holes, with opposite spins, are jumping back and forth between 
the central copper and its four surrounding oxygens. As a consequence of this there is an 
increased tendency to double occupancy in the central Cu. The reduction of the magnetic 
moment at the central Cu was also, found with qualitatively similar behavior, in Ref. P,p!3 



The coupling with the lattice takes advantage of a preexisting enhanced effective mass, 
and produces the self-trapping of a Zhang- Rice singlet. To follow the evolution of this singlet 
we define the following operator: Z] = (P-^-^d]^^ — P-^^dl^-^) where creates a hole with spin 
0" in a symmetric combination of the four oxygen orbitals. In Fig. |^ we show the expectation 
value of < Z'^Z > in the ground state as a function of A. This quantity is related to the 
probability of finding a ZR state in the cell. We can see a rapid increases of this value when 
the polaron is formed. However, the double occupancy also increases at the central Cu. 
The amount of double occupancy could be measured from its order parameter < nj^^nj ^ > 
computed at Cui. We found 25% of double occupancy when the polaron is formed, while in 
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the undeformed lattice this value is only 6%. From this, we conclude that the polaron is a 
ZR singlet but of a more general character that the one found in the undeformed lattice. 

V. SINGLE PARTICLE SPECTRAL FUNCTION 



In this Section, we analyze the evolution of the spectral weight for single-particle exci- 
tations when the coupling to the lattice is switched on. We first review the previous known 



results [p0| , pi| , p6| , p7| , for A = 0. There has been some theoretical and experimental contro- 
versy about how the Fermi level and the spectral weight evolve with doping in real materials 
and therefore we pay special attention to these issues, in the light of the present results. 



In particular, some experiments [E3] suggest that the chemical potential remains pinned in 



the gap of the insulator as the system is doped, whereas others [^], suggest that it shifts 
with doping as in an ordinary semiconductor. In both cases the gap tends to fill in with 
transferred spectral weight. 

The single-particle spectral function is defined as, 

9^A^) = 9l{^) + gfA^) (2) 

where 

9l (^) = E I < ^" + 1) 141^0 (iV) > r X 5 (a; - {N) + En{N + 1)) . (3) 

n 

To obtain gf^ [uj) we have to change and + 1 by Cjo- and — 1 respectively. 
\il)n (N) > is an eigenstate having eigenvalue En (N) where N is the number of holes; Eo{N) 
is the ground state energy. All the energies and gaps given below are in units of t. and g^ 
describe the photoemission (PES) and inverse photoemission (IPES) spectra respectively. 

Following Refs. |^ and , the expression for gi^j {uj) can be evaluated numerically 
using the Lanczos algorithm and the continued fraction expansion formalism 
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In Fig. H we show the single-particle spectral function {gi = gi^ + gii) at the Cu and O 
sites in the undoped (half-filling) case A = 0, i.e. for the undeformed lattice. The main 
features of this spectra have been identified in a previous study [^. First of all, we can 
distinguish the lower (LHB) and upper (UHB) Hubbard bands, which are well separated 
by Ud plus a typical hybridization energy. Below the UHB we find the nonbonding oxygen 
band (NOB) with a small weight on the Cu sites. The energy difference between these two 
bands is slightly larger than A + 2Upd- The quantity A + 2Upd is the (bare) charge-transfer 
gap expected from a ionic limit. 

Between the UHB and the NOB other structure emerges. We have indicated the mo- 
mentum of the states close to the Fermi energy. All of them have total spin 1/2. Due to 
the spatial symmetry of this lattice, the peak with momentum k = (0, 0) is absent but 



it reappears if tpp is switched on [^]. These highly correlated states (CS) correspond to the 



Zhang- Rice singlets band [^. The gap is given by the energy difference between UHB and 
CS. This energy difference is ~ 3 for our parameter set. 

When the system is doped with one hole and A = (Fig. |^) the Fermi level is shifted to 
/i+ = Eo(5) -Eo(4) |23|,g. This is exactly the position of the k = (tt, 0) peak in Fig. |. This 



can be contrasted with the HF mean field calculations []TT| where the stoichiometric system 
does not show spectral weight in the place where the chemical potential of the one hole doped 
systems sits. This is worth emphasizing because it shows that the HF results cannot be 
automatically invoked as an explanation of the photoemission experiments showing pinning 
of the chemical potential . The proximity of the next available state in Fig. |^ indicates 



that the system is metallic. The low energy spectral weight (LESW) between the Fermi 
level and the gap in Fig. ||, as defined in Ref. p3| , p4| , is equal to 1.85. This quantity should 
be equal to 1 for an ordinary semiconductor, and equal to 2 for the Hubbard model in the 



atomic limit. The value close to 2 is in agreement with the results of Ref. p3| , |2^ and it is 
due to strong covalency. It is worth noting that, at the HF level, departures from an integer 
value of LESW (in this case 2) are not allowed. Transitions from these states to scattering 



states above the Fermi level were identified |Tl| as the origin of the mid infrared absorption 
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in the optical conductivity. 

For < A < Ac, as discussed in the previous sections, the symmetry of the lattice is 
lowered and hence the charge tends to accumulate around the Cui and Cu4 sites (Fig. 1). 
This produces changes in the local g^s, but the total spectral weight remains similar to the 
undistorted one and the system stays metallic. 

As we have shown in Section III, the BLD becomes stable in the region A > Ac. In this 



case, the translational invariance is broken further. In Fig. |I0|, we show the spectral weight 
for A = 1.2. The self-trapping of the carrier around Cui makes the local hole-occupied g 
to increase. The same quantity is larger in Cu4 than in Cu2 indicating that there is still 
some amplitude for the CS state to hop out and in from the Cui, mainly to the next nearest 
neighbor's cell. The local g^s for the oxygens (not shown) around Cui (CU4) increase the 
hole-occupied (-unoccupied) spectral weight close to the Fermi level. Apart for this feature, 
they are similar to the total oxygen g [Fig. p!0|(e)]. 

The levels located slightly above the Fermi level in Fig. |T^ show a polaronic shift. This 
makes the fi~^ to sit, not in the position of the first peak in Fig. |^, like in the A = 
case, but somewhere in the gap. This is a direct consequence of the fact that we took an 
infinite mass for the O ions. If the mass were finite we expect Fig. ^ to be X- dependent. 
Accordingly, some intensity will show up at the position of the new /i+ corresponding to 
the polaronic configurations. However, for an infinite mass, such configurations can not be 
reached quantum mechanically by adding one particle to the half-filled case. Thus, their 
spectral weight vanish. It is tempting to associate this behavior with the experimental 
observation of growing spectral intensity and pinning of the chemical potential in the gap 
upon doping [^. The small spectral intensity in the undoped case due to a finite O mass 



will be hardly measurable and hence interpreted as the gap. Nevertheless the polaronic shift 
close to Ac is too small by itself to explain the effect. It can help to this interpretation when 
combined with other effects like impurities. The low energy spectral weight in this case is 
2.06. This increase is expected from a higher covalency but, in our case, contrary to Ref. 
23| , |2^ , this effect is a local one. For even larger values of A the polaronic shift increases 
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further in qualitative agreement with Ref. 
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VI. OPTICAL ABSORPTION 



In this Section, we report our results for the optical conductivity. The calculation of this 
quantity in the framework of exact diagonalization follows the procedure explained in detail 



in Refs. p0|,|29|] The optical absorption at {u) or the real part of the complex conductivity 



0" {(jj) is defined as: 

sites 

where 

The quantities: 

Jp = ir tjj+y ({m^}) ^c|^Cj+j,T-,(T ~ c|_|_yg.QT-,o-) (6) 

ia 

and 

are the current and kinetic energy in the direction of the applied field, respectively, r is the 
distance between two nearest-neighbors sites, in units of the lattice constant a. 

In Fig. |Tl|, we show the real part of the conductivity in the undoped case for A = 
(undeformed lattice). We can see a dominant feature formed by the peaks near u; = 5. These 
peaks correspond to a charge transfer of a hole from Cu to O ( UHB to NBO transition). 
As one would expect, this absorption feature appears at A + f/pd modified by a typical 
hybridization energy. 

As we have discussed in the previous section, the minimum single particle gap (between 
UHB and CS) is ~ 3. These transitions are forbidden in our cluster. 
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In the one-hole-doped system, and for A = 0, it is known that there is a shift of the Fermi 
energy towards the CS band. Consequently, we should see a peak at low energy. Indeed, 
this low energy structure may be observed in Fig. |T^(full line) where the first peak at ~ 2 
corresponds to transitions from CS to NBO. The structure around uj ^ 5 corresponds to 
UHB NBO transitions. The upper peak at ~ 9 corresponds to CS LHB transitions. 

Now, we examine the optical conductivity in the region where the BLD is stable. In 
Fig. |12| ( dotted line ), we show the optical conductivity in the case of A = 1.2. The most 
important change with respect to the undeformed case is a shift of the peaks to higher 
energies. Taking into account the previous analysis, we can associate this change to the 
polaronic shift of the single-particle levels. We also can see the set of peaks related to the 
UHB NBO transitions. These peaks are approximately in the same place as before. 
Another shift to higher energies can be seen in the peak associated with the CS LHB 
transitions. This peak appears at higher energies because it occurs between the self-trapped 
state and LHB. 



VII. SUMMARY AND CONCLUSIONS 



In this work, we have examined the self-trapping of an added hole in the Peierls- Hubbard 
Model using exact diagonalization in small clusters. In the presence of deformation, the hole 
forms a local singlet analogous to the Zhang-Rice singlet in the undeformed lattice but of a 
more general character, since double occupancy in Cu sites is not longer negligible. 

Our calculation is complementary to inhomogeneous HF studies and, in general, it 
confirms the qualitative behavior found in this approximation. However, there are important 
differences between the results obtained by these two approaches. Here we discuss these 
differences as well as the connection with other results in the literature. 

Firstly, there is a reduction of the magnetic moments of the Cu where the doped hole is 
self-trapped. However, this effect is much less dramatic than in the HF calculations. 
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Secondly, the chemical potential for zero A shifts with doping to a place where spectral 
weight already exists in the undoped compound p3| , p^ in contrast to the HF results. This 



is not the case for non-zero A. Such an anomalous behavior arises from the lack of lattice 
dynamics in both our calculation and in the HF approach. The observed polaronic shift does 
not seem strong enough by itself to explain the observed pinning of the chemical potential. 
We have shown that the spectral weight behaves in a similar manner as reported in pre- 



vious studies, where the coupling with the lattice was absent [^,3 . There is an anomalous 
transfer of spectral weight but, contrary to those studies, this transfer is related here to local 
effects, i.e. it manifests in local spectral functions. 

The peaks in the optical conductivity are also shifted consistently with the polaronic 
shift of the single-particle spectral functions. This behaviour is qualitatively similar to the 
HF plus RPA calculations [1^ at electronic energies in clusters of the same size. 

Thirdly, the small polaron appears in this study only above a critical value of electron- 
lattice coupling. This important difference with respect to HF calculations can be under- 
stood by the fact that these calculations neglect the translational motion of the correlated 
states and the associated kinetic energy. One expects that Ac will be determined by the 
competition between the kinetic energy and the polaronic effects. 

In this exact diagonalization study we have considered a quite small cluster. Unfortu- 
nately a determination of the size dependence is out of the present computational possibil- 
ities. However, if we assume that Ac scales with the bandwidth of the Zhang-Rice states, 
a strong size- dependence is expected since the effective bandwidth of strongly correlated 



systems, for small clusters, decreases with increasing size |27 



The physics close to Ac is highly nontrivial for a finite value of the ion mass The 
transition is no longer sharp but washed out by quantum fluctuations, and Ac becomes a 
crossover value. The renormalization of the effective mass of quasiparticles is moderate 
and not exponentially large like in the large A limit. In real materials, neither A nor Ac 
are known accurately. The fact that relatively lightly doped cuprates show already high 
mobilities coexisting with polaronic signatures, suggests that the system is close to Ac. A 
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related scenario is that the system forms a more extended polaron with a BED distortion 
in a hmited region, or fluctuates between both, localized and extended states PJT9|. 

In summary, we have given a complete characterization of the three-band Peierls- 
Hubbard model in small clusters in the adiabatic approximation using exact diagonalization. 
We have examined some properties which can not be properly considered within HF calcu- 
lations, which have been so far extensively used to study these systems. We believe that 
these effects can be generic to a wide class of materials. 
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FIGURES 

FIG. 1. Patterns of displacements of the oxygen ions along the Cu-0 bonds for: (a) Breathing 
Extended Deformation (BED) and (b) Breathing Local Deformation (BLD). 

FIG. 2. Ground state energy as a function of A for both, BED and BLD. The energy is 
measured with respect to the energy of the undistorted lattice. 

FIG. 3. Equilibrium lattice displacements of the O's along the Cu-Cu axis as a function of A 
for both type of deformations (BED and BLD). 

FIG. 4. Ac as a function of Ud- 

FIG. 5. Occupancies of the CUO4 cluster as a function of A in the BLD. We show occupancies 
in both Cui and Cu4 clusters. 

FIG. 6. < 5^ > as a function of A in the BLD. We show < S*^ > in Cui, Cu2 and CU4 clusters. 

FIG. 7. < Z^Z > asa function of A in the BLD. 

FIG. 8. Single particle spectral function for the half-filled case in (a) Cu and (b) O sites. 
The full and dashed lines correspond to PES and IPES experiments, respectively. The spectra is 
independent of A since for reasonable values of A the lattice is undistorted. 

FIG. 9. Single particle spectral function in (a) Cu and (b) O sites for the 1 hole doped case 
and A = 0. The full and dashed lines correspond to PES and IPES experiments, respectively. 

FIG. 10. Single particle spectral function for A = 1.2. The full and dashed lines correspond 
to PES and IPES experiments, respectively. We show the spectra in non equivalent Cu's (a,b,c) 
and the total spectra in Cu (d) and O (e). 

FIG. 11. Real part of the conductivity in the undoped case for A = 0. 
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FIG. 12. Real part of the conductivity in the one hole doped case for A = (full line) and for 
A = 1.2 (dotted line). 
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